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^ ■ ABSTRACT 

<D A new approach is proposed to phenomenological study of a generic unified 

supergravity model, which reduces to the minimal supersymmetric standard model. 
The model is effectively parametrized in terms of five low energy observables. In 
^ consequence, it is easy to investigate systematically the parameter space of the 

model, allowed by the requirement of radiative electroweak symmetry breaking and 
by the present experimental limits. 

Radiative corrections due to large Yukawa couplings and particle-sparticle mass 
splitting are included into the analysis and found to have important effects, in par- 
ticular on the degree of fine tuning in the model. 

In this framework there are presented the predictions of the model for various 
low energy physical observables and their dependence on the values of the top quark 
mass and tan/5 is discussed. Results are also given for the large tan/3 scenario, 
tan (3 ~ nit/nib. 

Our approach can be easily extended to non-minimal supergravity models, which 
do not assume the universality of the soft breaking parameters at the unification 
scale Mx- Such an extension will be particularly useful once the masses of some 
sparticles are known, allowing for a model independent study of the parameter space 
at M x . 
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1. Introduction 

The " supergravity inspired" minimal supersymmetric standard model (MSSM) 
is an interesting laboratory for studying various phenomenological aspects of super- 
symmetry. In this model the soft supersymmetry breaking parameters of the low 
energy lagrangian are given, via the renormalization group (RG) evolution, in terms 
of the unification scale parameters. These are: the universal gaugino and scalar 
masses M Q and m , respectively, and the trilinear scalar coupling A . Apart from 
the gauge and Yukawa couplings there are in the model two more free parameters: 
the Higgs superfield mixing mass \x Q in the superpotential and the corresponding soft 
breaking parameter B a . The actual mechanism of the soft supersymmetry breaking 
(and the origin of the //-term in the superpotential) are not yet understood theoret- 
ically. Although there are various theoretical suggestions with regard to the values 
of those parameters at the unification scale, they should be considered at present as 
free parameters of the model. 

One of the most attractive features of the supergravity models is the possibility 
of radiatively induced gauge symmetry breaking' 1 !, which becomes manifest at the 
tree level after the RG evolution of the parameters of the lagrangian from the unifi- 
cation scale M x to the scale M z . This way the hierarchy M z /M x becomes related 
to the presence in nature of large Yukawa couplings. The minimal model supple- 
mented by the requirement of radiative gauge symmetry breaking is a very predictive 
framework: with only 5 free mass parameters //<>, m , M Q , A and B Q (and the yet 
unknown top quark Yukawa coupling Y°), it interconnects 1) the 577(2) x U(l) 
gauge symmetry breaking with v =174 GeV, 2) the very constrained Higgs sector of 
the MSSM, 3) the sparticle mass spectrum and 4) the relic density of the lightest 
supersymmetric particle (LSP). In addition, it has recently been strongly reempha- 
sized that the predictivity of the model is greatly enhanced by the requirement of no 
fine tuning of the parameters^. Indeed, the requirement that the M WjZ be within 
the measured value without fine tuning of the unification scale parameters strongly 
limits the magnitude of the supersymmetry breaking scale and, in consequence, of 
the sparticle masses. 

In spite of the large number of papers devoted to the model, its analysis is not 
yet fully satisfactory (for a complete list of earlier references we refer the reader 
to the ref. [3]). Its full parameter space has not been systematically explored 
yet. Secondly, it has recently been understood that radiative corrections originating 
from the combined effect of large Yukawa couplings and particle-sparticle mass 
splitting have important effects on the mechanism of radiatively induced SU(2) x 
£7(1) symmetry breaking' 4 ! and on the Higgs sector' 5 ^. In consequence, one may 
expect important modifications in the predictions of the model. 

There are two main points in our strategy in this paper. Firstly, we set up the 
formalism to explore the full parameter space of the model. Our formalism is based 
on the bottom to top approach whose essence is as follows: we choose a set of values 
for 5 low energy physical parameters, and then look for the set of parameters at the 
unification scale M x which induces radiative gauge symmetry breaking and gives the 
chosen low energy values. Thus, we reverse the standard procedure which is to find 
the low energy parameters corresponding to a given set of the boundary conditions 
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at the unification scale. This way we can systematically study the phenomenology of 
the "supergravity inspired" MSSM compatible with the idea of radiatively induced 
symmetry breaking, with no additional assumptions (apart from their universality) 
about the values (and their correlation) of the parameters at Mx (for which we 
have little theoretical insight). In this approach it is straightforward to implement 
all the experimental constraints on the low energy parameters and to discuss sys- 
tematically the correlations between physical observables, which are predicted by 
the model. Effectively, the model is now parametrized by five low energy physical 
parameters. We would also like to stress that our approach is readily extendable to 
non-minimal supergravity models, which do not assume the universality of the soft 
supersymmetry breaking terms at the scale Mx, provided more than five low energy 
physical parameters is used as input. Once some sparticles are discovered and their 
masses known, such an extension will allow for a model independent study of the 
parameter space at the unification scale. 

Secondly, we consistently include radiative corrections generated by large Yukawa 
couplings and particle-sparticle mass splitting. To this end, we use the RG equa- 
tions for all parameters of the potential with heavy sparticles decoupled at their 
thresholds. In the context of radiative breaking of the 577(2) x U{1) symmetry this 
technique has already been advocated in ref . [6] . It has also been used to study simi- 
lar corrections to the Higgs sector in the MSSM and found to be quite accurate' 7 '. In 
this paper we use the RG equations of ref. [6] with squarks and gluinos decoupled at 
Msusy (the precise definition of this scale will be discussed in Section 2). Radiative 
corrections generated by the other sparticle-particle mass splitting are small and 
can be neglected. 

The question about the role of the large radiative corrections in the phenomeno- 
logical analysis of the model has already been given some attention in the literature 
(see ref. [3] and references therein). However, to our knowledge, no fully consistent 
treatment of such corrections at the 1-loop level, simultaneously in all four sectors 
mentioned above, exists yet. Their full impact on the phenomenology of the model 
has become easy to explore only within our bottom to top approach. 

The paper is organized as follows: In the next Section we present the details of 
our approach. Section 3 is devoted to a qualitative discussion of those predictions of 
the model which follow from the requirement of radiative gauge symmetry breaking. 
In Section 4 we discuss the problem of fine tuning of the parameters and emphasize 
the role of radiative corrections in sizably diminishing the degree of fine tuning. In 
Section 5 we present our results, after imposing the presently available experimental 
constraints. Finally, Section 6 contains a very brief discussion of the large tan/5 
scenario and a summary of this work is given in Section 7. 

2. The Bottom-Up Approach 

The details of our procedure are as follows. We begin with choosing, in addition to 
the "experimentally" known quantities 

a EM (M z ) = 1 /127.9 , sin 2 6 W (M Z ) = 0.2324 , 

m T (M z ) = 1.746 GeV, M z = 91.19 GeV , (2.1) 
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a set of values for the following low energy unknown parameters: 



M t , Ma, tan/5, uiq, 



(2.2) 



where Ma is the Higgs pseudoscalar mass and 771% and mfj are the soft (low energy) 
mass parameters of the squarks of the third generation. The physical squark masses 
are eigenvalues of the mass matrix. For the bottom to top approach this is the 
most economical choice of the low energy input parameters: Mz, tan/5 and Ma 
are directly related to the parameters of the tree level potential and the squark 
mass parameters determine the magnitude of the dominant radiative corrections to 
the scalar quartic couplings. So they must be fixed anyway at the begining of the 
analysis. To proceed further we have to discuss in some detail the unification of 
gauge and Yukawa couplings. 

2.1. Unification of gauge and Yukawa couplings 

The analysis of the unification of couplings in the MSSM have been addressed 
in a number of papers'^ , with the most recent and thorough discussion given in ref . 
[9,10]. At the 2-loop level, the RG equations for the gauge and Yukawa couplings 
are a coupled set of equations. However, in the calculation of the gauge couplings, it 
is a very good approximation to set the top quark Yukawa coupling to be constant ' 9 ' 
(for tan (3 ~ m t /mt„ the b-quark Yukawa coupling is as large as the t-quark coupling 
and the same procedure must be applied to both of them). Then the effect of the 
Yukawa couplings on the running of ctj is realized as a negative correction to the 
1-loop coefficients b^s (which is, anyway, small; see below). Effectively, for some 
chosen value of Y t , we can first discuss the gauge coupling unification. Starting 
from the experimental values for aEuiJ^z) and sin 2 9\y(Mz), we use the unification 
conditions to predict the values of the running strong coupling constant a 3 (Q) and 
the unification scale M x (for the purpose of this paper we neglect the experimental 
uncertainity in the determination of sin 2 0w)- After the inclusion of the 1-loop 
threshold corrections coming from the supersymmetric particles, we obtain for each 
gauge coupling the relation 



is the 1-loop threshold correction to l/ai(Q). The summation is over all sparticle 
and heavy Higgs doublet states with masses above Q; Vj is the contribution of each of 
those particles to the 1-loop beta function coefficient of the gauge coupling (when 
solving the RG equations at the 2-loop level, it is sufficient to consider threshold 
corrections at the 1-loop level). The Aj includes the threshold correction related to 
the top quark mass and the effect of the Y t on the running of c^. The importance 
of various terms in eq. (2.3) can be illustrated by considering their contribution to 




(2.3) 



where 7^ contains the 2-loop contributions to the beta functions, and 




(2.4) 
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the value of a^(Mz), obtained from the unification conditions. It is well known that 
the 2-loop corrections 7^ give ~10% increase in the value of a 3 (M z ) as compared 
to the 1-loop result and must be included into the analysis. The 2-loop top quark 
Yukawa coupling corrections to a 3 (M z ), included in A;'s, give —0.001, —0.0015, — 
0.002, -0.003 for F t 2 (M x )/47T ~ 0.1, 0.2, 0.4, 1.0, respectively! 10 ] and the top quark 
threshold correction is even smaller, (and of opposite sign) or at most similar, for 
M t < 200 GeV [9] . Altogether, for the range of the top quark masses and Yukawa 
couplings considered in this paper, those corrections do not exceed (1-2)% and we 
neglect them in our analysis. 

The contribution of the sparticle thresholds to l/a 3 (Mz) can be effectively 
parametrized in terms of only one parameter Tsusy as follows: 

fc(k)) = ^ 

where the scale Tsusy has the obvious interpretation in the (unrealistic) case of the 
fully degenerate sparticle spectrum: it is just the common mass of all the sparticles. 
Notice that, in this case, a non-zero threshold corrections to \/oLi{M z ) are present 
only if Tsusy > M z and, in consequence, the contribution (2.5) implies a decrease 
in the value of a 3 {Mz\ 

In the realistic case, in the presence of sparticle mass splitting, and in the ap- 
proximation in which one neglects the effects of the electroweak symmetry breaking 
on the sparticle masses, i.e. assuming no mixing, the scale Tsusy reads^ 10 ': 

Tsusy = m, ( f ( ^V'" ( (2.6) 

\ m 9 ) \ m d \ m Hj \ m Hj 

where nig, m g , mj, m^, and m H are the characteristic masses of the squarks, 
gluinos, sleptons, electroweak gauginos, Higgsinos and the heavy Higgs doublet, 
respectively. Here we have assumed that all the sparticle masses are larger than 
Mz (in the approximation of no mixing); if any sparticle mass is smaller then Mz, 
it should be replaced by M z when applying eq. (2.6). It is obvious that T S usy 
depends strongly only on the first two factors and in the model under investigation, 
with gaugino masses coming from a common supersymmetry breaking scale M Q at 
M x , Tsusy is well approximated by^ 10 ': 



w * ms Uw w "tej.' (27) 



Clearly, unless the Higgs mixing parameter in the superpotential \i > 600 GeV, the 
effective scale T S usy < M z and the sparticle threshold corrections tend to increase 
the value of a 3 . As it will be discussed later on, large values of /i are less natural and 
the most interesting range is 100 GeV < fi < 600 GeV. Therefore, for our model, 
the typical values are 15 GeV < Tsusy < Mz- The corresponding values of a 3 (Mz) 
obtaind from the unification condition vary from a 3 (Mz) = 0.123 for Tsusy = Mz 
(i.e. effectively, threshold corrections to a 3 (M z ) cancell out) to a 3 (M z ) = 0.127 
for Tsusy = 15 GeV (this case can be realized in practice e.g. by decoupling the 
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colour sparticles in the RG equations for cc^'s at m co i pa 3M Z and the remaining 
sparticles at M z (see eq, (2.6))). Experimentally, a 3 (M z ) is poorly determined and, 
conservatively, any value in the range (0.12±0.01) is still acceptable. We see from the 
above discussion that the sparticle spectra predicted by the model push us into the 
upper range of values for a 3 (Mz), unless we take less natural solutions. Of course, 
for full consistency, one should take care of the fact that the threshold corrections 
correlate the value of 0:3 with the specific solution obtained for the parameters of 
the model (\i in particular) which, in turn, slightly depend on the value of a 3 (M z ). 
However, there is an inherent uncertainity in the theoretical prediction for a 3 (M z ) 
which follows from the unknown contributions from the thresholds at the unification 
scale and from the higher dimension operators which are likely to be as important 
as the dependence on Tsusy^- Therefore, in this paper we take Tsusy = M z and 
use for the gauge couplings the supersymmetric RG equations down to the scale 
Mz- This gives a 3 (M z ) = 0.123, close to the central experimental value. For the 
unification scale we get M x = 2.4 x 10 16 GeV. 

Once the running of the gauge couplings is determined, we can proceed with 
the determination of the top, bottom and r Yukawa couplings in the whole range 
(Mz, M x )- In our procedure the Y T and Y t at M z are fixed by the values of the 
physical r mass M T and the chosen physical top quark mass M t and the value of 
tan/3. The running mass m t (M t ) is related to the physical mass (defined as the pole 
of the propagator) by the following 1-loop formula! 11 ] : 

m q(M q ) = I l (Mq) (2.8) 

For the running of m t (M t ) down to M z we use the 2-loop RG equations of the 

MSSM and we obtain Y t (M z ) from the equation m t (M z ) = Y t (M z )v/yfl + cot 2 /?. 
The running of the couplings Y t , Yt>, Y T in the range (Mz, Mx) is obtained by solving 
the 2-loop coupled RG equations with the boundary condition Y" b ° = Y° at M X - The 
solutions to the RG equations with those mixed (low and high energy) boundary 
conditions are obtained by iteration. For consistency with the RG equations for the 
gauge couplings we use the supersymmetric RG equations for the Yukawa couplings 
in the whole range (M z , M x )- Indeed, it has recently been shown' 10 ' that the same 
effective scale, Tsusy, governs the supersymmetric particle threshold corrections to 
the gauge and Yukawa couplings. 

The reason for taking the top and not the bottom mass as the input parameter is 
as follows: as soon as the top quark is discovered, its mass will be known with good 
precision whereas the physical bottom mass is not directly accessible to experiment 
and its value will remain very uncertain and model dependent. In this paper we 
present the results for the following sets of (M t , tan/3) values: a) 130, 2, b) 130, 
10, c) 180, 2, d) 180, 30, e) 130, 31 and f) 180, 51 (all masses given in GeV). It is 
interesting to comment on the predictions for the bottom quark mass obtained in 
each of those cases. They are summarized below: 
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m b (M z ) m b (M b ) M b 



b 



d 



c 



a 



c 



f 



4.01 5.57 6.08 

4.04 5.60 6.12 

3.30 4.71 5.16 

3.55 5.02 5.49 

3.81 5.33 5.82 

3.17 4.56 5.01 



The running mass m b (M b ) is obtained from m b {Mz) by means of the 577(3) x U(l) 
2-loop RG equations and the physical mass M b is given by the eq. (2.8). 

As said above, the physical bottom mass is a poorly known quantity but, still, 
following the particle data book it is expected to be in the range, say 



We see that the conditions of the gauge coupling unification and of the unification 
of the bottom and r Yukawa couplings result in too large a value of the M b , for 
moderately heavy top quark. It is so because the large value of a^iMz) (~ 0.12) 
gives too much of a strong correction to the running of m b (Mx) = m T (Mx), which 
has to be partly cancelled out by large, negative contribution from the top quark 
Yukawa coupling. In consequence, as discussed in detail in ref. [10], in this su- 
persymmetric grand unification scenario, the value of M b in the range (2.9) implies 
the top quark mass remarkably close to its quasi infrared fixed point value associ- 
ated with the "traviality" bounds on the top Yukawa coupling. A top quark mass 
M t « 130 — 140 GeV happens to be close to its infrared fixed point (i.e. consis- 
tent with (2.9)) only for values of tan/5 close to one, or very large values of tan/5, 
tan/5 > rrit(mt) / m b (mt) , above the value for which the unification of the three 
Yukawa couplings, Y° = Y h ° = Y°, occurs. Those values of tan (3 are the lower and 
upper (or even exeeding it) limits, respectively, for which radiative gauge symmetry 
breaking can be realized in the model and, at best, at the expense of a rapidly 
increasing fine-tuning of the parameters (to be discussed later on). The situation 
is different for M t ~ (180 — 190) GeV, which corresponds to its quasi infrared fixed 
point value for a large range of tan j3 values. 

In summary, exact unification of gauge and r and the bottom Yukawa couplings 
in the MSSM, together with the range (2.9) for the bottom mass, imply that M t < 
180 GeV can be realized only for tan/5 1 or tan/5 m t (m t ) / 'm b (m t ) , i.e. for values 
for which the mechanism of radiative gauge symmetry breaking can be realized only 
at the expense of a large amount of fine-tuning (or, in some cases not realized at 
all). From this point of view, the values M t > 180 GeV look more natural in the 
model. However, one should remember that the exact unification of couplings may 
be subject to the high energy threshold and higher dimension operator corrections, 
which can modify the above strong conclusions. This is true for both the gauge 
coupling and the r and bottom Yukawa coupling unification condition. For instance, 
if as(Mz) is smaller by A = —0.006 (the maximum magnitude of the high energy 
corrections, estimated in ref. [9]) then we get the following results: 



M b = 4.6 - 5.3 GeV 



(2.9) 
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m b {M z ) 


m b (M b ) 




a 


3.88 


5.30 


5.75 


D 


o.y i 




O.oU 


C 


2.97 


4.21 


4.60 


d 


3.42 


4.76 


5.18 


e 


3.70 


5.09 


5.53 


f 


3.08 


4.35 


4.74 



On the other hand, with a^iMz) = 0.123 and relaxing the bottom and r Yukawa 
coupling unification, for M b = 5.2 GeV one gets the following results for the ratio 
r = h b (M x )/h T (M x ): 





r 


a 


0.829 


b 


0.820 


c 


1.008 


d 


0.923 


e 


0.846 


f 


1.135 



The above examples illustrate the order of magnitude of the corrections to the exact 
unification conditions, which are necessary to make moderate values of the top quark 
mass and tan/? compatible with the range of values (2.9) for the physical bottom 
mass. 

In view of all the uncertainities discussed above, in this paper we present the 
results for all cases (a - f), with the gauge and Yukawa couplings fixed by the exact 
unification condition. Even in the presence in some of those cases of the likely 
inconsistency with the experimental value in the predicted bottom quark mass, this 
is a useful exercise since it allows to discuss interesting correlations present in the 
model, as a function of M t and tan/?. Moreover, modifications to a 3 (M z ) and/or 
to the relation Y T (Mx) = Y b (Mx), which are sufficient to bring down the value of 
M b into the range (2.9), have only little impact on the rest of our results. 

2.2. Determination of the parameters in the Higgs potential 

After those preliminaries we are now able to fix all the parameters in the Higgs 
potential at the scale M z 

V = mj | H l | 2 +m\ | H 2 | 2 -m\{e ab H*H b 2 + h.c.) 
+ 2^i I Hi | 4 +2^2 | H 2 | 4 +A 3 | Hi | 2 | H 2 | 2 
+A 4 | e ab El H\ | 2 (2.10) 

in terms of the values taken for the set (2.1) and (2.2). (As usual the boundary values 
of the parameters at Msusy read: m\ 2 = m 2 2 + /i 2 , m| = B/i where m 2 2 are the 
soft Higgs mass parameters, fi is the superhiggs mixing parameter and B is the 
corresponding soft breaking parameter, ei 2 = +1 and m| is defined to be positive). 
However, before we proceed further, it is useful to recall the conventional approach 
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to radiative SU(2) x £7(1) breaking. It makes use of the renormalization group 
improved tree-level potential evaluated at a renormalization scale Q ~ Mz- The 
running parameters in the potential, which are functions of the renormalization point 
Q, can be evaluated for any Q by solving the standard supersymmetric RG equations 
(written in the Landau gauge and in the (mass-independent) modified minimal 
subtraction scheme, MS) with some universal boundary conditions at the scale 
Mx- For Q ~ Mz, those RG equations sum up all the large logarithmic corrections 
proportional to ln(Mx/<3), which are reabsorbed in the running parameters of the 
potential (2.10). Due to the supersymmetric structure of the RG equations the 
tree-level potential maintains its supersymmetric form at any Q. Therefore, also at 
Mz, the parameters satisfy the supersymmetric relations: 

Ai = A 2 = i(^ + ^ 2 ), X 3 = \(9 2 2 - 9l), K = -\gl (2.11) 

This procedure neglects radiative corrections generated by the soft supersymmetry 
breaking (i.e. neglects logarithmic corrections of the type ln(m s to P /Q)), which are 
proportional to large Yukawa couplings and as we know now, can be very important 
and are of our concern here. One of the well known consequences is that the mini- 
mum of the tree-level potential is very sensitive to the choice of the renormalization 
scale Q. The choice Q = Mz gives, in general, unreliable results. The reason is that 
the loop correction AVi(M^) to the potential is not necessarily small.j] This has been 
demonstrated explicitly^: the sum V a (Q) + AVi(Q) is indeed scale-independent, up 
to 2-loop corrections which turn out to be small. For the sake of easy reference we 
give in the Appendix the 1-loop corrections to the parameters of (2.10): Amf(Q) 
and A\i(Q) (calculated in the symmetric phased i.e. before SU{2) x £7(1) break- 
ing) generated by the squark loops. This is the dominant part of the entire 1-loop 
contribution to the effective potential. 

For our purpose, the inconvenience of this approach is that one cannot work 
directly with the tree-level potential (2.10) and therefore the relations between the 
physical low energy parameters and the unification scale initial conditions are com- 
plicated. Fortunately, there exists the well-known equivalent method of including 
radiative corrections in question, which is based on the Appelquist-Carazzone de- 
coupling theorem. This is the tool to absorb in the running parameters of the 
tree-level potential (2.10) also the corrections generated by the particle-sparticle 
mass splitting. The main source of large corrections is in this case the top-stop 
splitting. At this point one should stress that sparticle thresholds contribute dif- 
ferently to different quantities and the dominant corrections to the scalar potential 
have different origin than the corrections to the gauge couplings. Therefore, to a 
very good approximation, the running parameters in (2.10) at the scale Mz can be 
related to their initial values at the unification scale by solving the RG equations in 
two steps: the supersymmetric equations in the range (Msusy, Mx) and the RG 
equations with squarks and gluinos decoupled^ at the common scale Msusy, in the 



3 The correct choice of the scale Q is such that AVi (Q) ~ but it is unknown a priori. 
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range (Mz, Msusy) ft By decoupling squarks at M S usy chosen to be equal to the 
low energy squark soft mass parameters one reabsorbes the logarithmic parts of the 
corrections (A. 7) to (A. 16) (in the approximation rag = mfj = m 2 D ) and, of course, 
also higher order leading logarithmic corrections, into the running parameters of the 
potential (2.10). As long as we work in the region (Mz, Msusy) with the 1-loop 
RG equations, this is a consistent approximation: all non-logarithmic corrections in 
eqs. (A.7-A.16) are formally the effects of the same order as the neglected 2-loop 
effects in the RG equations. However, earlier results' 4 ' strongly indicate that the 
former are much more important than the latter. In particular, as discussed earlier, 
they are important for the potential to be renormalization scale independent (in 
our case, this means, independent of the choice of Msusy in a reasonable range 
around the soft squark masses). Also, it is known that the 2-loop effects in the RG 
equations for the evolution of the quartic couplings are very small' 12 '. Therefore, 
in this paper we include also the non-logarithmic part of the 1-loop corrections 
(A.7-A.16) by properly choosing the matching conditions for the parameters m 2 
and Aj at Msusy- More specifically, we fix the scale Msusy so that the contribution 
(A. 8a) to Am 2 vanishes. All 1-loop corrections to the quartic couplings are then 
included by replacing at M S usy the continuous matching conditions given by the 
relations (2.11) by the conditions A, — > A, + AAj with AAj's given by the eqs. (A. 13- 
A.16) (everywhere in this paper we use the approximation m 2 ^ 2 = m 2 ^ 2 = m? Di = 
m 2 D = 2mg 3 — mfjj which follows from the supersymmetric RG equations after 
neglecting terms proportional to the electroweak gauge couplings). Similarly for 
m 2 , i — 1, 2, 3 , we take discontinuous matching conditions at Msusy '■ rnf + Am 2 
where Am| is given by the eq. (A. 9), Am 2 . - by the (A. 8b) and Am 2 is given by the 
expression (A. 7b), always with A = M^^y-I The technical way of implementing 
the discontinuous matching conditions for m 2 's is discussed after eq. (2.17). 

After the decoupling of the sparticles the largest 1-loop corrections to the tree 
level potential obtained after the evolution are given by the top quark loops. The 
top quark should therefore be decoupled at A top chosen to minimize those 1-loop 
corrections to the vevs Vf 



4 The mass difference between the third and the remaining two generations of squarks can be 
neglected as the latter contribute little, due to the small Yukawa couplings. Gluino mass is, in 
most cases, close to the squark masses. The other sparticles are kept in the spectrum down to 
Mz- The reason for using in this paper the RG equations with only squarks and gluinos decoupled 
is that, on the one side, radiative corrections to the scalar potential from the other sparticles are 
very small and, on the other side, we want to analyze the light sparticle spectrum predicted by 
the model. However, as far the scalar potential is concerned, very similar results are obtained by 
using below Msusy the RG equations of the standard model with two Higgs doublets. 

5 We neglect the part (A. 7a) of Am\ which is very small for our choice of Msusy- One should 
also remember that always m± 3> m\ and the correction Amf is relatively much less important 
than Am|. 



dAVj 

dv { 







where 
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We get then A top ~ M t /y/e and for the interesting range of the top quark masses 
we can use the approximation A top w M Z - Thus, by the RG evolution in the 
range (Mz, Msusy) one reabsorbes in the tree level parameters the remaining large 
corrections to the potential. 

2.3 Bottom— up approach 

The procedures outlined above are used at two different stages of our approach. 
First, we come back to the problem of fixing all the parameters in the Higgs potential 
(2.10) at the scale Mz in terms of the values taken for the set (2.1) and (2.2) and use 
this procedure to obtain the couplings \(Mz) from their boundary conditions (2.11) 
at Msusy. In general \(M Z ) = \{M U M T , M su sy, tan/3, g k , (M x )) as their RG 
evolution depends on the Yukawa couplings. Given the values of \{Mz) we obtain 
the mj(Mz), % = 1,2,3 , in terms of Mz, tan/3 and Ma(Mz) by minimizing the 
potential (2.10)0. The relevant formulae can be found eg. in ref. [13] and, for easy 
reference, are collected in the Appendix A. 

The second stage of our approach consists in finding the soft supersymmetry 
breaking parameters m , M , A , B Q and the parameter /x c at the Mx which give 
radiatively induced spontaneous gauge symmetry breaking with the chosen set of 
values for the low energy observables (2.1) and (2.2). To this end, we notice that 
in both regions (M z , M SUS y) [20] and (M S usy, M x ) [6] the RG equations for the soft 
Higgs and gaugino mass parameters as well as for the trilinear couplings A^s can 
be schematically written as follows (t = In -rpr): 



(The equation for rh\ will be discussed later. Note that the parameter B appears 
only in the equation for rh\.). The coefficients C lk in eqs.(2.12) depend on the gauge 
and Yukawa couplings and change at M SUSY (below M SUSY they also depend on 
the Aj couplings). Defining the multidimensional vector 



(where eg. AiAj denote all possible products of A^s present in the model and 
similarly for the other quantities) and using eqs. (2.12) one has 



6 We neglect the running of the pseudoscalar mass Ma, i.e. we work in the approximation 
Ma(Mz) = Ma(Ma) which is correct to at least few GeV accuracy. 




Cmm m i + ^mM + ^rnA + ^ k — 1,2, Q,U .. . 



(2.12) 



X = {MiMj, M t Aj, AiAj, ml /i 2 ) 



d 
dt 



X{t) = C(t) X{t) 



(2.13) 
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(where C(t) is a linear operator) i.e. we have a set of linear, homogeneous equations. 
Its solution reads: 

X(t) = R(t,t ) X (Q 
with R (t,t a ) = T exp^C (s) ds^ (2.14) 

where T means chronological ordering in case of non-commuting operators C(t). In 
view of the boundary conditions at Mx'- 

A° = A , M° = M , m° = m 

the solution (2.14) can be written more explicitly as[](we recall that rhf = mf + 
/<'-• i = l,2). 

m\(M z ) = Cij(M t , Mr, gl tan/3, M SUSY ) P°, i = 1,2 (2.15) 

where 

P°={Ml ml & Al M A ). 

Quite similarly one gets 

m 2 QtU (M SUSY ) = V Q>v (M t , M T , gl tan/?, M SUSY ) P° (2.16) 

The coefficients c^- and 6q v in the algebraic equations (2.15) and (2.16) are 
found by integrating numerically the RG equations with the boundary conditions 
like P° = 1, P° = for i / 1 etc. This is a two-step procedure: first, the 
super symmetric equations are integrated from Mx to Msusy with those boundary 
conditions at M X - Then the RG equations with decoupled squarks and gluinos are 
integrated from Msusy to Mz, with the boundary values at Msusy obtained from 
the first integration. Msusy is fixed in advance by the choice of the squark masses, 
eq. (2.2). 

Thus, we have a set of 4 algebraic equations for 4 independent parameters 
M Q , m , Ho, A Q . Finally, B Q is fixed from 

ml(MsusY) = KM SUSY ) B(Msusy) 

by the equation 

ml(MsusY) = d x [L M + d 2 [i A + d 3 fx B (2.17) 

(as before, the coefficients di are calculated from the RG equations). 

At this point we can complete our discussion of implementing the matching con- 
ditions for Am- at Msusy- All the terms (A. 7b), (A. 8b) and (A. 9) are proportional 
to the parameters \i and A^ taken at M S usy- It follows from the eqs. (2.12) that, 
analogously to (2.15), 

H = Cfl 

A{ = A Q + c M M 

7 Analogous equations have been discussed earlier in ref. [14] with the coefficients derived ana- 
lytically (in the case of no sparticle decoupling and negligible Yb). 
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and, therefore, the corrections Am- have identical structure in terms of the unifica- 
tion scale parameters as the mf's themselves. So, their inclusion amounts to some 
modification of the coefficients q,- and bq^ in the eq. (2.15) and (2.16). 

Our procedure is now completed and can be summarized as follows: for any 
chosen set of values for the physical parameters (2.2) we find the corresponding low 
energy parameters of the effective lagrangian and then solve the algebraic equations 
(2.15), (2.16) and (2.17) to find the lagrangian parameters at M x which give those 
chosen values of tan/3, M A , m 2 Q and m 2 v . We note that for each set (2.2) there are 
up to 8 solutions for the parameters at Mx- the sign of m is taken to be positive by 
convention (the lagrangian is symmetric under an overall change of sign of all soft 
breaking terms) but the solutions appear in pairs of ±M G and ±/x (the solutions 
e.g. for M D < are the same as the solutions for M Q > and the reversed signs 
of all parameters, with exception of m ). The remaining two- fold ambiguity follows 
from the fact that the set (2.15-2.16) is equivalent to a quadratic equation for the 
ratio A /M . This overall eight-fold ambiguity is simply a reflection of the fact that 
the set (2.2) is not a complete set of observables in the model: different solutions 
give different prediction for other observables. Finally, one should stress that for 
a given set of values for the observables (2.2) our equations may have no solution 
at all. This is, in fact, the most often encountered case for a random choice of the 
values (2.2) and it simply means that the given set cannot be obtained in the model. 
Our procedure in this paper is to fix M t and tan /3 and then to scan systematically 
the values of Ma, win and mfj in the range, say, from to 1 TeV, by varying them 
in steps of a few GeV in search for solutions to our equations. Thus, for a given 
choice of M t and tan/3 we get a range of values for Ma, rn 2 q and mfj which are 
obtainable in the model (i.e. for which the solutions to eqs. (2.15) and (2.16) exist). 
For each solution we can then calculate the values of the other physical observables. 
Finally, we impose on the obtained solutions the "no fine tuning" constraint and all 
the presently available experimental constraints. 

3. Qualitative Features of the Solutions 

In this section we discuss those properties of our solutions which follow from 
the requirement of radiatively induced gauge symmetry breaking alone. We also 
include the constraint of all masses squared to be positive and the usual one of no 
deeper lying minima of the scalar potential that break charge and/or colour[15]: 
A 2 < 3(mQ + m 2 , + ml) and similarly for the remaining Aj's at low energy, as well as 
the analogous condition at the unification scale. Finally, we demand that the scalar 
potential should be bounded from below at the scale Mx- 

The " naturalness" constraint and the constraints from the present experimental 
limits will be discussed in the next sections. As an introduction to our results let 
us recall certain properties of the low energy potential (2.10). One can verify that 
radiative corrections stabilize the potential l 6,14 !. So, the low energy mass parameters 
rhi, i=l,2,3, which correspond to the chosen values of tan/3 and Ma, do not need, 
in general, to satisfy at Mz the well known stability condition for the supersym- 
metric tree level potential: m\ < {m\ + m 2 ) 2 /4:. Secondly, as already mentioned 
earlier, up to the standard model 1-loop corrections, the effective potential is now 
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renormalization scale independent in the range (Mz, Msusy)- This means that 
the gauge symmetry is broken at the tree level already at M S usy, i-e. at the scale 
where the quartic couplings in the potential (2.10) recover their supersymmetric 
values.^ Note, however, that the individual parameters in the potential and in the 
lagrangian, in general, run sizably in the range Mz, Msusy (particularly the quartic 
scalar couplings !). Physical quantities at the scale Mz receive in addition the SM 1- 
loop corrections (particularly, from the top quark and vector boson exchange) which 
are included into the analysis by our choice of the renormalization scale Q = Mz- 
Next, it is instructive to analyze the coefficients in eqs.(2.15) and (2.16) and to 
understand them qualitatively with help of the RG equations. A sample of coeffi- 
cients for mf, rriq and mfj is given in Table 1. Of special interest is their dependence 
on the top quark mass, tan (3 and their modification by radiative corrections. A re- 
mark is in order here: it is the usual procedure followed by many authors to solve 
the supersymmetric RG analytically in some approximation. We could, of course, 
follow the same approach here to calculate the coefficients in eq. (2.15) and (2.16). 
However, the approximations involved in this procedure make it very difficult, if not 
impossible, to discuss the most interesting, just mentioned, dependences. Therefore, 
we prefer to base our discussion on numerical results, guided by general qualitative 
features which can be inferred directly from the RG equations. Let us first look 
at m\ and m\ f° r sma U t an P- Taking into account the definitions rrif = mf + /i 2 
and the initial conditions mf = m 2 , we expect from the RG equations (dominated 
by the terms proportional to Y t - see the Appendix B) that for m\ the coefficients 
are close to 1 for m 2 and /x 2 , positive and smaller for M 2 and negligible for A 2 a and 
M A . For m\ the coefficient of /z 2 is expected again to be 0(1) but the coefficient 
of m 2 gets additional negative contribution from the terms proportional to Y t which 
largely cancels the +1 contribution from the initial condition. The most interesting 
is the coefficient of M 2 which is 0(1) and negative ! This looks contradictory to the 
positive gauge terms in the RG equation for m\ (the same as for m 2 ) and is the 
reflection of the running of the scalar masses in the equation for m\. We have 

y 9\Ml + ... 

0, M 3 (M X ) = M (3.1) 

and this strong dependence on M 2 is transmitted to m\ through the negative term 
—3Yt(mQ + mfj + . . .) in the equation for m\. Naively a second order effect, it is in 
fact the dominant one due to the large logarithmic enhancements (and the running 
of Ms) and it plays an important role in understanding our results. In particular 
it explains strong dependence of the M 2 coefficient in m\ on Msusy'- the larger 
is this scale the smaller is the absolute value of the negative M 2 coefficient, due 
to the earlier decoupling of the gluino term in eq.(3.1). This behaviour has also 
important implications for the fine tuning problem (to be discussed in the next 
section) before and after radiative corrections. Indeed, in our formalism, neglecting 

8 Another way of seeing this is to remember that the driving force for the difference in the mf 
and m\ RG evolution from the scale Mx (which gives the breaking at low energy) are squarks and 
they are decoupled at Msusy- 
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radiative corrections to the scalar potential means setting Msusy — Mz for any 
values of the stop mass parameters and we see (the numbers in brackets in Table 1) 
that the M 2 coefficients in m\ are then, for the same set of values for the variables 
(2.2), much larger in the absolute values than after radiative corrections. 

It is also very important to recall that the low energy values of ttlq and mfj are 
almost directly given by the values of M%. The next-to-leading dependence is on 
w? but with the coefficients by factor O(10) smaller then for Ml. Again, the M 2 
coefficients decrease substantially after radiative corrections. 

The tan/3 and M t dependence of the coefficients is mainly determined by the 
behaviour of the Yukawa couplings and their contribution to the RG equations. For 
fixed M t and increasing tan f3 the Yj, becomes larger and the y t -smaller (eventually, 
Y t /Y b = 1 for tan (3 = m t (M t ) / 'mf,(M t )) . Accordingly, the coefficients in the equation 
for m\ behave as follows: the Ml coefficient is driven from positive to negative 
values, the one for the m 2 Q is also decreasing but remains positive, due to the ml 
contribution to m\ , and the M A coefficient is slowly rising from zero to ~ 0.4. This 
behaviour follows from the increase of Y b which starts to play the same role for m\ 
as Y t for m\. At the same time the Y t is decreasing (but it changes much less than 
Y b ) and, in consequence, the coefficients in the equation for ml change just in the 
opposite direction. Still, even for tan (5 ~ m t /m b the coefficients remain different 
enough from those in the expression for m\ to assure the existence of solutions ' 16 1. 
The dependence on M t (for fixed tan/3) can be easily understood in a similar way, 
and follows from changing Y t . (We note also the decrease of the nl coefficients with 
increasing Y t ). 

Generically, the mass parameters of the low energy potential satisfy the relation 
m\ ^> m\ and m\ < 0. Thus, to understand qualitatively the pattern of our 
solutions, in the first approximation we can take m\ = 0. In a large region of the 
parameter space we then have: 



with b, c, a, 7 positive and a and (3 generally smaller and of both signs. This 
structure, together with the relation Ma ~ m\ and 



explains the correlations seen in our solutions for the absolute values of the pairs 
(M G , fi ), (M a , m Q ) and (/i , m Q ). It is most visible in the (M a , /i ) correlation 
(or, equivalently, (M 2 , //)) which is interesting in the context of the neutralino 
properties predicted by the model. It is clear from the coefficients in Table 1 that, 
as long as we can neglect the term M A in the equation for m\ (and neglecting also 
the other small coefficients), crudely speaking, the solutions exhibit approximately 
linear correlation Ml (7/a)/^. The actual allowed regions in the space (M 2 ,/i) 
are determined by the integration over the soft scalar masses mq and mfj, and 



over the pseudoscalar mass Ma- They are shown in Fig. 5 for tan/3 = 2 and 10 
for M t = 130 GeV and for tan/5 = 2 and 30 for M t = 180 GeV. A remark 



m\^ a Ml + b m 2 + c fi 2 
ml ~ -aMl + (3m 2 + 7/^ w 



(3.2) 




(3.3) 
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has to be made here. In this section we refer to all the solutions shown in the 
Figures. The distinction between experimentally excluded or allowed and natural 
or unnatural solutions (which is marked in the Figures) will be discussed in the 
subsequent sections. The dependence of the (M 2)/ u) correlation on tan/5 and M t 
follows directly from the discussed above properties of the coefficients and in the first, 
crude, approximation is just the reflection of the changing slope a in M 2 m afi 2 ,. The 
correlation (M 2 ,/i) has direct consequence for the predicted neutralino properties. 
It is well known^ 17 !, and follows directly from the structure of the gaugino-higgsino 
mass matrix, that in the limit of large M 2 3> fi the lightest neutralino becomes 
almost pure higgsino with M n = and in the limit of large fi 3> M 2 - pure 
gaugino with M n « M 2 /2. The pattern observed in the (M 2 ,/i) plane explains 
qualitatively the results shown in Fig. 7. We plot there the gaugino content Z n = 
Z 2 (l, 1) + Z 2 (l,2) of the lightest neutralino as a function of the neutralino mass, 
as predicted by our solutions. (As usual we define xi — l)B + Z(1,2)W 3 + 
Z(l,3)Hi + Z(1,4)I?2 where Z(i,j) is the real orthogonal matrix that diagonalizes 
the neutralino mass matrix.) In particular, the M t and tan/3 dependence of those 
results is easily interpretable in terms of the analogous dependence of the (M 2 ,/x) 
correlation. 

The above simple discussion breaks down and a new feature of our solutions 
becomes visible for the values of M t and tan /3 such that the coefficient of the M 2 
term in the equation for m\ is much smaller in its absolute value than that of the 
M A term. For instance, this is the case for M t = 160 GeV and tan (3 = 20. Then 
the dependence on M Q of the rh\ comes from the M A term and is correlated with 
the values of A a . Positive A Q allows for larger positive M a , since then the M A term 
in m\ (no longer negligible) is positive and helps to cancel out the negative M 2 term. 
Neutralinos with large gaugino content (large fi, small M 2 ) appear mainly for light 
ttiq (hence small M Q and (— \A\ ) branch) and large Ma « rh\ (hence large /i). For 
even larger values of M t and tan/3 (e.g. M t = 180 GeV, tan/3 = 30) the coefficient 
of the M 2 term in the equation for m\ becomes again large in the absolute value 
(but negative) and the effect disappears. 

Another interesting correlation present in the model is seen in the (mj , M Q ) plane 
{mi is the mass of the lightest squark), as illustrated in Fig. 10. It follows from 
the structure of eq. (3.3) and from the qualitative features of the coefficients shown 
in Table 1: the mj is almost directly related to Mj^\. One can also expect some 
correlation in the (mf, Ma) plane. As explained above, M a is correlated with \x a 
(up to the contribution of the m and A Q terms) and therefore M\ = m\ should be 
correlated with mj. However, as seen from Fig. 12, the dependence on the additional 
parameters makes the final correlation rather weak. An obvious correlation exists 
between m T and the m , Fig. 11. The qualitative difference between Fig. 11A, 
C and Fig. 11 B, D is due to the negative contribution of the Y& term to the RG 
equation for m?, which is increasing for large tan/3. 

In Fig. 1-4 we also show our solutions in the planes (m , M a ), (m , A Q ), (A Q , M ) 
and (A /m , B /m ). Clear positive correlation is observed in the (m , A Q ) projec- 

9 Strinctly speaking, this is true for the soft squark masses. However, the results obtained after 
diagonalization of the stop mass matrix (presented in Fig. 10), with the mixing terms predicted 
by our solutions, show very similar pattern. 
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tion. Interesting structure is seen in the (A /m , B /m ) plane: there is an almost 
linear correlation (see eq. (2.17) and Table 1) and two branches of solutions exist, 
corresponding to the two solution with /i = ±|/x|, and with the gap depending on the 
value of m\. Only weak correlation is observed in the (m , M ) plane. Obviously, it 
will be very interesting to study correlations in 3 and more variables. 

So far we have been discussing the qualitative features of the predictions based 
on the solutions to eqs.(2.15), (2.16) and (2.17), i.e. following from the requirement 
of radiatively induced gauge symmetry breaking. In the next two sections we impose 
two additional constraints: no fine tuning and the present experimental limits. 



4. Fine Tuning 

The basic feature of the model is to predict the low energy parameters in terms 
of the unification scale parameters which, hopefully, can be obtained in future from 
some underlying theory. This procedure involves extrapolations of many orders 
of magnitude in the energy scale. In this context it is important to address the 
question of the degree of fine tuning of the parameters at Mx which is necessary to 
get the correct low energy results. Precise values of the input parameters do not, 
by themselves, mean fine tuning yet. Only when the precision we need for the input 
parameters is much higher than the precision we require for the low energy results, 
it means fine tuning. In other words it is the question of response in the low energy 
results Qi to small changes in the high energy input Pj, which is measured by^l 

Aij J^lf^Il (4.1) 

Equivalently this is also a measure of cancellations in equations like (2.15): for 
instance, calculating the mass parameters rhf in the low energy lagrangian in terms 
of the high energy parameters Pi (using the coefficients like those in Table 1) the 
single terms on the rhs should not be much larger than the calculated masses. To 
some extent it is a matter of taste how big cancellations (i.e. how big A^-'s) are 
still tolerable. The coefficients Cy in eqs. (2.15)— (2.17) depend on large logaritms 
In M x /M z = 0(30), so we take here the attitude that cancellations of this order of 
magnitude are still "natural". 

Our actual interest is to avoid fine tuning when calculating the physical quantities 
like Mz (or v), tan/5, M t , etc. So we look for solutions with acceptable all the 
following derivatives 

Pi dQ 

Ay = -7^- , Qi = {Y t ,Y T ,v, tan (3, M A , m Q ,m v } 

P % = {Y t °, Y T °, M , m , n , A , B } (4.2) 

We shall consider a solution to be "natural" if 

| | < A , for i,j = 1,...7 (4.3) 

where, following the arguments given above, we sort the solutions according to 
A = 10, 30 and no constraint on |Ay|. It turns out, however, that in the range of 
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the low energy parameters studied in this paper, Ma, tuq, mjj < 1 TeV almost 
all solutions satisfy |Ay| < 100 (for moderate values of tan/5 discussed here). 

For given values of v and tan/3 the constraints (4.3) mean upper limits on the 
unification scale parameters, or equivalently, upper limits on Ma and the stop masses 
for which "natural" solutions exists. In general, it is difficult to avoid large cancel- 
lations simultaneously in all the low energy mass parameters. 

There are two main effects which ease the problem of fine tuning after radiative 
corrections. These are: increase of the absolute value of ml for the same set of 
low energy parameters (see eq. (A. 4) and remember that, generically m\ < 0) and 
(the discussed earlier) decrease of some coefficients in eqs.(2.15)-(2.17), particularly 
the ones of the Ml terms. The overall effect is that the maximal acceptable Mi 
considered as a function of the max |Ay| is, with radiative corrections included, 
almost twice as large as without. Similar effect exists for m- t and /i. This results in 
a sizable enlarging of the "natural" regions (for any chosen value of A) on all the 
Figures. For the illustration of this point we mark in Figs, 1, 5, 7, 12 (by the solid 
contours) the "natural" regions, for A = 30, obtained without radiative corrections. 
One should also remark that for A = 10 the allowed parameter space of the model 
is small even with radiative corrections included. The fine tuning problem becomes 
more severe for tan/5 close to 1 and to m t /m b , i.e. close to the limits of the tan/? 
range obtainable in this model. Solutions for the large tan (3 case will be discussed 
in more detail in Section 6. 

A comment is in order here on our method of calculating derivatives OQi/dPj 
with Pk fixed, for k ^ j. Note that the coefficients c^- in eqs.(2.15)-(2.16) are fixed 
for given values of the low energy parameters (in particular they depend on the 
decoupling point which is fixed in terms of the low energy squark mass parameters). 
Thus, it is straightforward for us to calculate the inverse matrix dPj/dQi,Qk fixed 
for, k i, and to find by inversion the matrix of interest. The direct calculation 
of the latter is in our approach cumbersome because fixing the unification scale pa- 
rameters rather then the low energy parameters makes the coefficient c^- dependent 
on them (the decoupling point depends then on Pi's). Thus, we calculate: 



dQ l 




' dP k 




dPj 




dQi 





Note, that to get the correct answer for the lhs derivatives one has to calculate the 
complete matrix on the rhs. Most often, the largest derivatives are those of the 
v with respect to Y t , fi Q and A Q and of tan/3 with respect to B D , m , A D and Y t . 
Finally, one should remark that in eq. (4.2) we consider only the parameters which 
are used in the calculation as independent input parameters. It also happens that 
the derivatives of some other physical quantities (e.g. sparticle masses), taken with 
respect to the high energy set P i: become large, even when Ay's defined in eq. (4.2) 
are small. 
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5. Experimental Constraints and Final Results 

Let us recapitulate our procedure: For some chosen values of the low energy phys- 
ical parameters M t , tan/5, Ma, mg and my (see eq. (2.2)) we look for solutions 
to eqs. (2.15)- (2.17). These are the equations for the unification scale parameters, 
derived from the requirement of radiatively induced gauge symmetry breaking in 
the minimal supergravity model. The existence (non-existence) of solutions for a 
given set of input values means that this particular set can (cannot) be realized in 
the model. In this paper we present results for a few generic top mass and tan/3 
values: M t = 130 GeV and tan/3 = 2 and 10; M t = 180 GeV and tan/3 = 2 
and 30. In the next Section we discuss very large tan/3 cases: tan/3 ~ m t /mb. In 
each case the Ma, itiq and my have been scanned in the region from to 1 TeV 
in steps of a few GeV. Once some solutions for the unification scale parameters 
are found we can predict the values of any other low energy observables, corre- 
sponding in the model to our set of input data. We present our results in form of 
the scatter plots obtained by our scanning procedure. In this paper we limit our- 
selves to two-dimensional scatter plots, but multi-dimensional correlations are also 
of interest. In Fig. 1-4 we present our solutions in the unification scale variables 
(m G , M Q ), (m G , A ), (A Q , M a ) and (B /m , A /m ). Projections onto the vari- 
ables (M 2 , n), (M n , Ma), (Z n , M n ), (M c , M A ), (M c , M n ), (m h M Q ), (m f , m ), 
(m t ~, Ma) and (M c , m t ~) are shown in Fig. 5-13. (We recall the notation: M Q , m Q , fj, Q , 
A Q , B Q are the high scale parameters; M 2 is the low energy SU (2) gaugino mass, fx is 
the low energy supersymmetric Higgs mixing parameter; Ma, M n , M c , m~ t , m f are 
the pseudoscalar, and the lightest neutralino, chargino, stop and stau physical 
masses, respectively, and Z n is the gaugino content of the lightest neutralino. Three 
different classes of solutions are introduced: a) solutions already excluded by the 
present experimental constraints, b) allowed by the present experiments and sat- 
isfying the criterion max|Ajj| < 30, c) allowed by experiment and with 30 < 
max |Ajj | < 100. (In the considered region of Ma, rriQ, my < 1 TeV and for those 
moderate values of tan/3, solutions with max|Ajj| > 100 occur quite rarely.) It 
is worth pointing out that, in our 2-dimensional projections of the 5-dimensional 
parameter space, those three different classes of solutions are often placed in over- 
lapping regions. For clarity of the Figures, only part of the obtained solutions is 
plotted: only as many of them as is needed to properly mark various domains in 
the parameter space. Additional important information given in some of the Fig- 
ures is the sensitivity of our results to radiative corrections. The solid contours in 
Fig. 1,5,7,12 encircle the regions where the class (b) solutions are present, with no 
radiative corrections included. In the same Figures, the dashed contours show the 
regions where max |Ajj| < 10. 

For the purpose of this paper we have taken the following "experimental" limits 
as the constraints for our solutions: 

m- g > 120 GeV, m t > 45 GeV, M c > 45 GeV 

\Z 2 (1,3) -Z 2 (l,4)| <0.03, \Z(1,3) Z(2,3) -Z(l,4) Z(2,4)| < 0.04 
(for M n < 45 GeV) (for M ni + M n , 2 < 45 GeV) 

mg > 100 GeV (for the first two generations of squarks), 
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m~ t > 45 GeV, M A > 40 GeV. 



(5.1) 



Some of the above limits are somewhat rough representation of the actual experi- 
mental data. In particular, strictly speaking, at present there is no absolute lower 
limit on the stop mass' 18 !. The bounds from pp colliders do not apply to it; those 
bounds assume degenerate squarks and this assumption fails for strong L R mixing 
of the t-squarks. Also searches at the e + e~ colliders are affected since the L-R mix- 
ing can reduce the Zt\t\ coupling^. In the model considered, each solution gives 
definite predictions for mass mixings and couplings and in principle a more detailed 
phenomenological analysis with respect to the constraints from the present pp and 
e + e~ data is possible for each solution. We postpone it, however, for future study 
and take (5.1) as sufficient for our discussion here. 

As it has already been mentioned several times, the model shows a high degree of 
correlation among various physical quantities which results from the small number 
of free parameters. We have discussed in Section 3 that a number of correlations 
follows from the requirement of radiatively induced gauge symmetry breaking. We 
can add now to this discussion some remarks following from the fine tuning and 
the experimental constraints. As expected, the larger the values of the high en- 
ergy parameters the more fine tuning is required (see Fig. 1-4). For instance, for 
max \Aij\ < 30 we get M < (200-400) GeV and m < (400-900) GeV, depending 
on the values of the top quark mass and tan /3. Simultaneously, the ratio m /M can 
vary, approximately, in the range (0.1-10). As we know from Section 3, the gauge 
symmetry cannot be broken for m rs M q rs 0. At least one of them must be visi- 
bly different from zero. The experimental constraints put somewhat stronger lower 
bound on M Q (> 70 GeV) than on m D (> (20 — 30) GeV) and acceptable solutions 
with m « A Q w B Q « exist. 

At the quantitative level, the predictions of the model show strong dependence 
on the top quark mass, tan (3 and on the inclusion of the radiative corrections. For 
a detailed survey of the predictions for the physical (low energy) observables we 
refer the reader to Fig. 5-12. Here we would like to stress again the interesting 
correlations between M 2 and //, between M t and Z n (the gaugino content of the 
lightest neutralino), tan/3 and Z n , Z n and M n , mi and Ma etc. Similarly to the re- 
quirement of radiative gauge symmetry breaking, the experimental constraints (5.1) 
have "global" effects, i.e. each individual constraint affects not only this concrete 
quantity but also all the others through the rigid structure of the model. An in- 
teresting example is the interplay of the experimental lower bounds on the gluino, 
chargino and neutralino masses. In the model the three masses are interconnected 
by the parameters M 2 (or M ), \i and tan/3 and it varies from case to case which 
experimental bound is most operative in fixing e.g. the effective lower bound for the 
neutralino mass. For neutralinos with large higgsino content (/i small, M 2 large), 
i.e. with large coupling to the Z, the LEP data give directly M n > 45 GeV. For 
gaugino-like neutralinos (M 2 small, /i large) the bound on M n « M 2 /2 follows from 
the bound on the gluino mass (M 2 Ri 0.3 nig > 40 GeV) or from the LEP bound on 
the chargino mass (for M 2 <C // and fi > My/ one has M c ~ M 2 ). More precise num- 
bers follow from the well known complete chargino and neutralino mass matrices 

10 Similar remark applies also to stau leptons. 
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and depend on the actual values of the parameters in the given case. As seen e.g. 
in Fig. 7, the lower bound for M n is indeed in the range (20-50) GeV, depending 
on M t and tan (3. 

Another interesting example is the effective lower limit for the pseudoscalar 
mass Ma- We have taken Ma > 40 GeV as the experimental lower limit for 
Ma but, as seen e.g. in Fig. 6,8 and 12, the effective bound (which follows from 
the other constraints (5.1)) is much higher, Ma > 175, 125, 205, 150 GeV for 
(M t , tan 0) = (130, 2), (130, 10), (180, 2) and (180, 30), respectively. One can 
understand it as follows: M\ ~ m\ and, for the discussed here moderate values of 
tan/3, all coefficients in the equation (2.15) for rh\ are positive and large (see Table 
1). Therefore, small M\ implies small values of M%, m 2 and /x 2 ,. However, as 
discussed above, there are lower bounds on those parameters from the experimental 
lower bounds on the gluino, neutralino and chargino masses which, as shown in Fig. 
6 and 8, give the effective lower bound for Ma of the order O(200 GeV). 

For larger values of M t and tan/5 (e.g. M t = 180 GeV and tan/5 = 30) the 
coefficients in m 2 of the M 2 and \j? Q terms become very close to each other and of 
opposite signs. The above argument does not then work any more: arbitrarily large 
values of M 2 ps /i 2 can give small Ma- However, we see that, in the case considered, 
the coefficients of the M 2 and /x 2 terms in the m\ are very different and, since m\ ks 0, 
to obtain radiatively induced gauge symmetry breaking we need // 2 ^> M 2 (as seen 
in Fig. 5). Hence, the terms M 2 and /i 2 in the m\ cannot cancel each other and 
we again get the effective lower limit on Ma higher then the present experimental 
lower limit. Also note that similar arguments explain the correlation between light 
neutralinos, with large gaugino content, and heavy pseudoscalars which is observed 
for M t = 180 GeV. 

For even larger values of tan/5 the situation again looks differently: as will be 
discussed in more detail in the next section, for tan/5 ~ mt/mj, all the coefficients 
in the expression for m\ are very close to the coefficients for m\ and a light pseu- 
doscalar A becomes compatible both with radiative symmetry breaking and with 
the experimental constraints. 

In general, we recall that the relative magnitude of the coefficients of the M 2 and 
// 2 terms in the expressions for m\ and m\ is responsible for most of the structure 
which follows from the mechanism of the symmetry breaking alone. The lower limit 
on Ma is an example of its interplay with experimental constraints, for changing 
values of M t and tan /5. 

In Table 2 we present several examples of the generic sparticle spectra obtained 
in this model. In Table 3 we also give results obtained with and without inclusion 
of radiative corrections for the same sets of the input parameters. We see that the 
difference between the two sets can be very important. 

6. Solutions for Large tan (3 

It has been discussed for some timet 21 ' that the observed pattern of the quark 
masses, namely strong isotopic spin symmetry breaking in the third generation elec- 
troweak doublet, can be understood as the property of the vacuum. In models with 
two (or more Higgs doublets) such that two different vevs v 1 and v 2 are driving 
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the down - and up - quark masses, it is conceivable that tan/3 = v 2 /vi = m t jm^ 
with the top and bottom Yukawa couplings (almost) equal. Such a scenario is par- 
ticularly interesting in the context of supersymmetric grand unified models (they 
require at least two Higgs doublets in the low energy effective lagrangian) which, 
if based on the groups SO{\$$) or E 6 with simple Higgs structure, predict Y t = Y\, 
at the unification scale. This scenario has recently been strongly advocated in a 
detailed study of the fermion mass and mixing textures^ . 

An interesting question is whether the minimum of the potential with large 
tan/? can be obtained by radiative gauge symmetry breaking. The answer to this 
question is indeed positive, as shown for the first time in ref. [16] and subsequently 
investigated in a number of papers^. 

Our bottom to top approach is particularly suitable for a systematic search for 
large tan (3 solutions and a study of their low energy properties. A sample of such 
solutions is shown in Fig. 14-17 for M t = 130 GeV, tan/3 = 30 and 31, and 
M t = 180 GeV, tan/? = 50 and 51. The most important points about the large 
tan (3 solutions are the following: The degree of fine tuning is rapidly increasing when 
we approach the value of tan (3 such that Y° = For instance, for M t = 130 GeV 
and tan/3 = 31 (the corresponding value of the Y t °/Y£ at M x is 1.04) we have not 
found any solutions with max |Ay| < 100 whereas for tan = 30 (Y t °/Y£ « 1.08) there 
exist solutions with max |Ay| < 30. This is easily understandable from inspection of 
the RG equations: the masses m\ and m\ run now very similarly and in consequence 
the coefficients for m\ and m\ in eqs. (2.15) are very close to each other. As an 
example, for M t = 130 GeV, tan (3 = 31 and M SU sy = 460 GeV we have: 

m \ = -0.95 M 2 - 0.27 A 2 + 0.41 m 2 + l. 20 /i 2 + 0.47 M A 
m j = -1.07 M 2 - 0.10 A 2 + 0.45 m 2 + 1.20 /i 2 + 0.39 M A 
ml = 0.18 nM a - 0.41 n A + 1.08 /i B 

At the same time, the input values of m 2 's at Mz are: m\ ~ (6—9) x 10 3 GeV 2 , m\ ~ 
—5 x 10 3 GeV 2 and m\ « 70 — 160 GeV 2 . Clearly, we need large values of the 
parameters at M x (in particular of M )[]] and a sizable amount of fine tuning to 
solve those equations. 

The above example illustrates several other points typical for the large tan/3 
scenario. These are: m\ ~ (see eq. (A. 5) and small pseudoscalar mass M\ m 
ml + m\. The acceptable solutions exist now for small values of Ma, since for large 
tan p all the coefficients in the equations for m\ and m\ are close to each other and 
moreover, in the combination M\ = m 2 + m\ the terms M 2 and /z 2 tend to cancel 
each other. Therefore, the arguments of the previous Section do not apply: large 
tan (3 solutions in general give light pseudoscalar Higgs bosons and a much heavier 
spectrum for squarks and gluinos. 

Another point is that m\ ~ can be realized in two different ways: small \l or 
cancellations between the three terms. It is the latter possibility which is chosen by 
most of the solutions. 

11 The small differences in the coefficients are driven by the residual difference in the Yukawa 
coupling and also by the different U(l) charge assignement for the right-handed up- and down- 
squarks. The latter (and the slepton contributions) are sufficient for the existence of solutions even 
for Y° = Y£ but at the expense of very large values of M and max |Ay | ~ O(10 3 ). 
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Finally, we stress that large tan (3 solutions are particularly sensitive to radiative 
corrections included in our calculation. 

In conclusion, large tan f3 scenario is compatible with radiative gauge symmetry 
breaking, although the exact equality Y° = can be achieved only at the expense 
of very strong fine tuning. On the other hand, allowing for, say, (5-10)% difference 
in the two Yukawa couplings at Mx brings the necessary degree of fine tuning into 
acceptable range. It is quite possible that such difference is generated by high energy 
thresholds and contributions from higher dimension operators, but their impact on 
the Yukawa coupling unification has not been investigated so far. 

7. Conclusions 

In this paper we have proposed a bottom-up approach to unified supergravity 
models which give the minimal supersymmetric standard model as the effective low 
energy theory. In this approach the supergravity models are parametrized in terms 
of a number of low energy physical parameters. The minimal model, with universal 
boundary conditions at the unification scale and five free mass parameters, has been 
studied in detail. With more than five low energy input parameters, our set of 
equations can be used to investigate non-minimal models, which do not assume the 
universal boundary conditions at Mx- This will be particularly interesting once 
some sparticles are discovered and their masses known. 

Even in the minimal model, there are several virtues of the bottom-up approach. 
One is that the number of free parameters is reduced to four (plus the top quark 
mass) since the mass Mz is known. Secondly, a systematic and complete investiga- 
tion of the low energy predictions of the model can be easily performed, by scanning 
the space of four unknown low energy input parameters. Finally, with a suitable 
choice of the input parameters, our approach is most convenient for inclusion of ra- 
diative corrections generated by large Yukawa couplings and particle-sparticle mass 
splitting: by choosing the soft stop masses as the input parameters the dominant 
1-loop radiative corrections are precisely known directly in terms of the input pa- 
rameters. 

The complete inclusion of radiative corrections generated by large Yukawa cou- 
plings and particle-sparticle mass splitting is the other main point of this work. 
They are found to have strong effects on the degree of fine tuning needed in the 
model and on the final results. 

In the above framework, low energy phenomenology of the minimal model has 
been explored in some detail. In particular, we have imposed all the presently 
available experimental constraints on the predictions of the model and obtained 
complete information on those regions of the low energy physical observables which 
are interesting for future experimental exploration. 

Due to the very rigid structure of the minimal model, both the requirement of 
radiative gauge symmetry breaking and the experimental constraints have strong 
"global" effects: they introduce a lot of correlations between different quantities 
predicted by the model, as discussed in the previous sections. Particularly interesting 
and illuminating is the dependence of the predictions on the values of the top quark 
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mass and tan/3. In this respect, our result based on numerical solutions to the RG 
equations show a very rich structure, which is difficult to see fully in appproximate 
analytic calculations. On the list of interesting results are, among others, neutralino 
as the lightest neutral sparticle, the gaugino content of the lightest neutralino, the 
lower limit for the pseudoscalar mass Ma and, in general, correlations among masses 
(and couplings) of sparticles. The model predicts a rich sparticle spectrum to be 
accesible at the next generation of colliders. 

In the subsequent paper' 23 ' we calculate in the same approach the relic abundance 
of the lightest neutralino and impose the additional constraint VLh 2 < 1 on the 
predictions of the model. 

As a further extension of the present work, it will also be very interesting to 
explore systematically other constraints which may follow from more specific as- 
sumptions about physics at the grand unification scale, such as e.g. proton life time 
constraints in the supersymmetric SU(5) model' 24 '. 
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APPENDIX A 

The minimization of the potential (2.10) gives 

tan 2 (3(rh\ + \ 2 v 2 ) = m\ + X 2 v 2 + (A x - X 2 )vf 



2mi 



ml + ml + Aif \ + \ 2 v 2 + (A 3 + A 4 ) f 2 sin 2/3 



(A.l) 



where tan/5 = v 2 /vi and i> 2 = v \ + t> |. 
The pseudoscalar mass M A reads: 



777.1 



7 2 + rn 2 + Aiw 2 + A 2 ^2 + (A 3 + A 4 )v 2 

« m 2 + 7772 + AA 2 ^2 ( A - 2 ) 

The last relation holds in the approximation in which we neglect all but the dominant 
radiative correction, from the top-stop mass splitting, to A 2 . For 777 2 's we get the 
following relations in terms of the input parameters M A , tan/3 and v 2 : 

, tan 2 (5 — 1 



777? = Ml sin. 2 (3 + \iv 2 - - 2 



tan 2 + 1 



- (Ai + A 3 + A 4 )tj 2 



M - Sm/3+ 2 M| ia^TT 



(A.3) 



7772 = M\ COS 2 /3 - X 2 V 2 



tan 2 /3 



1 + tan 2 
M 2 A cos 2 /3 - -M 2 z tan H~] - AX 2 v 



- cos 2 P(X 3 + X 4 )v' 



2 z tan 2 /3+l 



, tan 2 /? 
1 + tan 2 p 



(A.4) 



m 2 = -Ml sin 2/5 



(A.5) 



1-loop corrections from the squark loops to the parameters of the potential (2.10) 
(calculated from Feynman diagrams in the "symmetric" phase of the model ' 6 ') are 
as follows: 



m l( e ff) = ™i + Am h 



1,2,3 



where 



(47T) 2 Am 2 = 3 £ \\gld A + Yf (a (m 2 QA ) + a (m 2 DA 



+ H 2 Yff(m 



2 2 



+ A 2 DA Yff(m 2 QA ,m 2 DA 



(A.6) 

(A.7a) 
(A.7b) 



(An) 2 Amj = 3 £ f" I^a + (a D (m^) + a G (m 2 ^)) 



(A.8a) 
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+ V 2 Yff Kv<J + A uX 2 f « A ,m> UA )} (A.8b) 



Z A 

+ A DA Yff(ml A ,ml A )) (A.9) 



where 



/ (mlml) = -1 + 1 log -| + 2 2 2 log -f A.10 
^ 1 ^ to? - to| A 2 to| - to? A 2 



d A = 2a (to 2 ^) - a (to^ a ) - a (rojj (A.ll) 

and 



/ 2\ ^ \i f d d k i 2 / m 2 

2 



(27r) d A; 2 - m 2 



r/= — - + 7£ -log47r (A.12) 



We recall the notation: m 8 = mf + /j 2 , « = 1,2, where m 2 's are the soft super- 
symmetry breaking Higgs boson masses and \i is the supersymmetric Higgs mixing 
parameter; to| = Bfi; Qa are left- handed squark doublets (the index A numbers 
the generations), Ua and Da are the up and down right-handed squark SU(2) sin- 
glets, respectively; the trilinear soft supersymmetry breaking scalar couplings are: 
Vl = Yi A A La , yj) = Y^A Da and = Y^A Ua where F's denote respective Yukawa 
couplings; g^s are gauge couplings. 

The scale A = M S usy has been chosen so that the term (A. 8a) vanishes (see the 
text). 

1-loop corrections to the quartic Higgs couplings are as follows: 
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4Ai = g\ + g% + 



2(4tt) 



E 



8 4i 



- [ -9t- 8 9 2 1 Yf+2AYf 



log 



m 



A 2 



+ (3 g\ + \gt + 24 F/ 4 - 12 fe 2 r/ - 4 g\ Y f) log 



m, 



Qa 
A 2 



(A.13) 



4A 2 



<?2 + 



2(4tt) 2 



E 



+ (|^-16^ a +24^ 4 ) log 



A 2 



'3^ + ^ + 24^ 



log 



m 



Qa 



A 2 



(A.14) 




2A 4 = 



(An)' 



E 



3 4 
2^2 



3^ 2 (^ 2 +F/))log 



m; 



Qa 

A 2 



+ 6y/y/(iog^ + /( 



^A'^A, 



(A.16) 
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APPENDIX B 



For easy reference we list here the 1-loop supersymmetric renormalization group 
equations for the Yukawa couplings and for the soft supersymmetry breaking parame- 
ters! 19 ]. 

2 jY* = Y* fa + 3gj - 3Y^ - £ (if 2 + 3Y/)) (B.l) 

2 Jt Y * = Y * (f 91 + \ dl + 392 ~ 3Yf ~ Y * ~ 5 + 3r/2 )) (R2) 



2 j t Y U A = Y U A ( f g! + H g l + 3 , 2 - 3 y/ - Y/ - 3 £ ^ ) (B.3) 



| m| A = 3^ 2 2 M| + ^ 2 M 2 - Y^ 2 (m| A + m\ A + m? + A 2 J (B.4) 
| m| A = AglMl - TYf (m\ A + m\ A + m\ + A 2 J (B.5) 



d 9 16 

■ m; 



' u £ 3 2 M 3 2 + 3g 2 M 2 + 1 g\M\ - Yf [m% A + m 2 DA + m\ + A^) 



-Y/ (m 2 ^ + + m\ + A 2 ^) (B.6) 



| ™L = f flM + \ g\Ml - TYf (m\ A + tu 2 Da + m\ + A 2 Da ) (B.7) 
| = y ^ + y £ 2 M 2 - 2Y/ (m^ + m 2 UA + m\ + A 2 ^) (B.8) 



-m 2 = 3^M 2 + ^ 2 M 1 2 -3 5:Y/(m 2 3A +m^+m 2 + ^ 

-£^ 2 R. + ™L + ™ 2 + ^L) (b.9) 



| m 2 = 3 y 2 M| + ^M 2 - 3 £ ^ K + m u A +m * + A U (B.10) 
^ = ^ [^1 + ^ - E fr* + 3Yf + 3Y* 2 )) (B.ll) 
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±B = - (3glM 2 + gf Ml ) - £ (A Lb Y° 2 + 3A Ds F/ + SA^Yf) (B-12) 

| A LA = - (3, 2 2 M 2 + g 2 M^J - ZA\ A Yf - £ (^if 2 + 3^ s F/ ) (B.13) 

= - (f sgM 3 + 3<7 2 2 M 2 + I^M X ) - A^y/ - 3A DA Yf 

-T.i^Y^ + SA^Yf) (B.14) 



It 



16 13 

y g 2 3 M 3 + 3g 2 2 M 2 + — g\M, 



- 3A UA Yf - A Da Y 



A 2 



(B.15) 



d / M, 



= 



d 
~dt 



9i = Wi 



(B.16) 



(B.17) 



Here t = hi -^f and 6 X = 1 + 10n G /3, 6 2 = 5 — 2n G , 6 3 = — 9 + 2n G . The above 
equations are used in the range (M X) M SUSY ). In the range (M SUSY , M z ) (i.e. 
after decoupling of squarks and gluinos we use the equations given in ref. [6]. 
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TABLE 1 



Coefficients in Eqs. (2.15-2.17) for M SUSY = 300 GeV. 
The numbers in brackets are without radiative corrections (M S usy = M z ). 
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TABLE 2 



Generic spectra of the supersymmetric particles with the lightest 
pseudoscalar A or neutralino or chargino (underlined). 
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Table 3 



Comparison of the supersymmetric particle spectra, 
for the same sets of the input parameters, 
with and without (the numbers in brackets) radiative corrections 
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FIGURE CAPTIONS 



Fig. 1. The regions in the unification scale parameters (m , M Q ), obtained by de- 
manding radiative gauge symmetry breaking, in form of the scatter plots ob- 
tained by the scanning procedure described in text. These are 2-dimensional 
projections of the solutions in the 5-dimensional parameter space. Crosses: 
solutions already excluded by experimental constraints discussed in the text; 
squares: solutions compatible with experimental constraints with max |Ay| < 
30; dots: solutions compatible with experimental constaints with 30 < max | A y - 
< 100. Dashed contours: solutions with max |A y -| < 10. Solid contours: solu- 
tions compatible with experimental constraints and with max |Ay| < 30, but 
with no radiative corrections. 

A) M t = 130 GeV, tan (5 = 2; B) M t = 130 GeV, tan/? = 10; 
C) M t = 180 GeV, tan/3 = 2; D) M t = 180 GeV, tan/3 = 30. 
In case C) there are no solutions corresponding to the contour criteria. 

Fig. 2. Same as Fig. 1 for (m , A Q ). 

Fig. 3. Same as Fig. 1 for (A M ). 

Fig. 4. Same as Fig. 1 for (B /m , A /m ). 

Fig. 5. Same as Fig. 1 for the low energy parameters (M 2 , /i). 

Fig. 6. Same as Fig. 1 for the physical lightest neutralino mass M n and the pseu- 
doscalar mass M A . 

Fig. 7. Same as Fig. 1 for the gaugino content of the lightest neutralino, 
Z n = Z 2 (l, 1) + Z 2 (l, 2) and its mass M n . 

Fig. 8. Same as Fig. 1 for the lighter chargino mass M c and the pseudoscalar mass 
M A . 

Fig. 9. Same as Fig. 1 for the lighter chargino mass M c and the lightest neutralino 
mass M n . 

Fig. 10. Same as Fig. 1 for the physical lightest stop mass mi and the universal gaugino 
mass at M x , M . 

Fig. 11. Same as Fig. 1 for the physical stau mass m f and the universal scalar mass 
at Mx, rn . 

Fig. 12. Same as Fig. 1 for the stop mass mi and the pseudoscalar mass Ma- 

Fig. 13. Same as Fig. 1 for the chargino mass M c and the stop mass m- t . 

Fig. 14. Solutions for M t = 130 GeV and tan/3 = 30. Notation as in Fig. 1 (dots 
correspond now to solutions with 30 < max|Ay| < 250). The corresponding 
value of Y t /Y b at M x is 1.08. 
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Fig. 15. Solutions for M t = 130 GeV and tan/3 = 31. There only exist solutions with 
max | Ay | > 30. Plotted are solutions with max|Ay| < 500. The correspond- 
ing values of Y t /Y b at M x is 1.04. 

Fig. 16. Same as Fig. 15 for M t = 180 GeV and tan/5 = 50. The corresponding values 
of Y t /Y b at M x is 1.18. 

Fig. 17. Same as Fig. 15 for M t = 180 GeV and tan/5 = 51. The corresponding values 
of Y t /Y b at M x is 1.12. Plotted are solutions with max |Ay| < 750. 
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